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Abstract. In this note, we take up the study of weak convergence for stochastic dif- 
ferential equations driven by a (Liouvillc) fractional Brownian motion B with Hurst 
parameter H € (1/3,1/2), initiated in [3J. In the current paper, we approximate the 
d-dimcnsional fBm by the convolution of a rcscalcd random walk with Liouvillc's ker- 
nel. We then show that the corresponding differential equation converges in law to a 
fractional SDE driven by B. 



1. Introduction 



The current article can be seen as a companion paper to [3], to which we refer for a 
further introduction. Indeed, in the latter reference, the following equation on the interval 
[0, 1] was considered (the generalization to [0, T] being a matter of trivial considerations): 

dy t = a(y t )dB t + b(y t )dt, y = aeR n J (1) 

where a : ~R n — > ~R nxd , b : M™ — ► W 1 are two bounded and smooth enough functions, and 
B stands for a <i-dimensional fBm with Hurst parameter H > 1/3. 

Let us be more specific about the driving process for equation ([T]): we consider in the 
sequel the so-called <i-dimensional Liouville fBm B, with Hurst parameter H £ (1/3, 1/2). 
Namely, B can be written as B = (B 1 , . . . , B d ), where the B l, s are d independent centered 
£T) ' Gaussian processes of the form 

B\ = \ {t-r) H --*dWl (2) 



O 



O ■ Jo 

o ; 

Q\ ! for a c?-dimensional Wiener process W = (W 1 , . . . , W ). This process is very close to the 

usual fBm, in the sense that they only differ by a finite variation process (as pointed out 
in p]), and we shall see that its simple expression (j2J) simplifies some of the computations 
in the sequel. In any case, B falls into the scope of application of the rough paths theory, 
which means that equation (CQ) can be solved thanks to the semi-pathwise techniques 
contained in [6j [7J [10] . The natural question raised in [3] was then the following: is it 
possible to approximate equations like ([1]) in law by ordinary differential equations, thanks 
to a Wong-Zakai type approximation (see [HI H31 E] for further references on the topic)? 

Some positive answer to this question had already been given in [5] , where some Gauss- 
ian sequences approximations were considered in a general context. In [3], we focused on 
a natural and easily implementable (non Gaussian) scheme for B, based on Kac-Stroock's 
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approximation to white noise (see [H1II2])- However, another very natural way to approx- 
imate B relies on Donsker's type scheme (see [TT] for the case H > 1/2 and [I] for the 
Brownian case), involving a rescaled random walk. We have thus decided to investigate 
weak approximations to (CQ) based on this process. 

More precisely, as an approximating sequence of B, we shall choose (X £ ) e>0 , where X 6 ' 1 
is defined as follows for i = 1, . . . , d: consider a family of independent random variables 
{Vk'i k > 1, 1 < i < d}, satisfying the 

Hypothesis 1.1. The random variables {r]l; k > 1, 1 < % < d} are independent and 
share the same law as another random variable 77. Furthermore, rj is assumed to satisfy 
E (77) = 0, E {rj 2 ) = 1 and is almost surely bounded by a constant k n . 

We then define X e,t in the following way: 

X i ' £ {t)= [ {t + e-r) H ~^e l ' £ {r)dr, (3) 
Jo 

where 

-, +00 

^):=-E^MA)(^ (4) 

k=l 

Notice that X £ is really a process given by the convolution of the rescaled random walk 
9 e with Liouville's kernel. 

Let us then consider the process y e solution to equation (CQ) driven by X £ , namely: 

dyZ = *(y*)dXZ + b(yZ)dt, y £ = aeM n , te[0,T]. (5) 
Our main result is as follows: 

Theorem 1.2. Let (y e ) e> o be the family of processes defined by and let 1/3 < 7 < H, 
where H is the Hurst parameter of B. Then, as e — > 0, y £ converges in law to the process 
y obtained as the solution to (QJ|, where the convergence takes place in the Holder space 
CT([0,l];E n ). 

This theorem is obtained invoking many of the techniques introduced in [3] . In the end, 
as explained at Section [21 most of the technical differences between the two articles arise 
in the way to evaluate the moments of quantities like f(r)9 h£ (r) dr for a given Holder 
function /, and to compare them with the moments of Gaussian random variable. This is 
where we shall concentrate our efforts in the remainder of the paper, mostly at Section [3j 

2. Reduction of the problem 

We shall recall here briefly some preliminary steps contained in [3] , which allow to reduce 
our problem to the evaluation of the moments of a specific type of Wiener integrals. 

First of all, we need to recall the definition of some Holder spaces, in which our con- 
vergences take place. We call for instance Cj([0, 1]; M d ) the space of continuous functions 
from [0, lp to Mr, which will mainly be considered for j = 1 or 2 variables. The Holder 
norms on those spaces are then defined in the following way: for / G C2QO, 1]; K d ) let 

||/|| M = sup J£*L and ^([0,l];M rf ) = {/eC 2 ([0,l];M d ); < 00} . 

s,te[o,T] |£ — s r 

The usual Holder spaces ([0, 1]; M. d ) are then determined by setting \\g\\^ = \\8g\\n for a 
continuous function g G Ci([0, 1]; IR d ), where Sg G C2QO, 1]; K d ) is defined by 5g st = gt~ g s - 



We then say that g G C± ([0, 1]; ~R d ) iff \\g\\n is finite. Note that || • || M is only a semi-norm 
on Ci([0, 1]; M. d ), but we will work in general on spaces of the type 

Cf >a ([0, 1]; R d ) = {g : [0, T] -> V; ^ = a, |MU < oo} , (6) 

for a given a G V, on which is a norm. 

The second crucial point one has to recall is the natural definition of a Levy area for 
Liouville's fBm: 

Proposition 2.1. Let B be a d- dimensional Liouville fBm, and suppose that its Hurst 
parameter satisfies H G (1/3, 1/2). Then 

(1) B is almost surely a ^-Holder path for any 1/3 < 7 < H. 

(2) A Levy area based on B can be defined by setting 

dB u ® dB v , i.e. Bl(i,j)= dBl dBl, i, j G {1, . . . , d}, 

J s J s J s 

for < s < t < T . Here, the stochastic integrals are defined as Wiener-Ito integrals when 
% 7^ j, while, when i = j , they are simply given by 

J*dBi£dBi = ~(Bi-Bi) 2 . 

(3) The process B 2 is almost surely an element of C| 7 ([0, 1]; M dx ' i ), and satisfies the 
algebraic relation 

B 2 - B 2 SU - Bl = (B u - B s ) ® (B t - B u ) , 
for allO < s < u < t < 1. 

These algebraic and analytic properties of the fBm path allow to invoke the rough path 
machinery (see [HI H CD]) in order to solve equation ([T]): 

Theorem 2.2. Let B be a Liouville fBm with Hurst parameter 1/3 < H < 1/2, and 
a : M n — > M. nxd be a C 2 function, which is bounded together with its derivatives. Then 

(1) Equation admits a unique solution y G C^W 1 ) for any 1/3 < 7 < H, with the 
additional structure of weakly controlled process introduced in [7j . 

(2) The mapping (a, B, B 2 ) ^ y is continuous from W 1 x Cj(R d ) x C^(M. dxd ) to Cj(R n ). 

One of the nice aspects of rough paths theory is precisely the second point in Theo- 
rem [221 which allows to reduce immediately our weak convergence result for equation (CQ), 
namely Theorem II .2\ to the following result on the approximation of (B,B 2 ): 

Theorem 2.3. Recall that the random variables rf k satisfy Hypothesis and let X e be 
defined by For any e > 0, let X 2,£ = (X^f(i, j)) s ,t>o;i,j=i,...,d be the natural Levy's 
area associated to X s , given by 

ft 



X%(i,j)= / (X^-X^)dX^, (7) 

J s 

where the integral is understood in the usual Lebesgue-Stieltjes sense. Then, as e — > 0, 

(X £ ,X 2 - £ ) ^ (B,B 2 ), (8) 

where B 2 denotes the Levy area defined in Proposition ^. 1\ and where the convergence in 
law holds in the spaces Cf(M d ) x C^iM.^^, for any // < H. 



The remainder of our work is thus devoted to the proof of Theorem 12.31 

As usual in the context of weak convergence of stochastic processes, we divide the proof 
into weak convergence for finite-dimensional distributions and a tightness type result. 
Furthermore, the tightness result in our case is easily deduced from the analogous result 
in®: 

Proposition 2.4. The sequence (X e ,X 2 ' e ) £>0 defined at Theorem \2.3\ is tight in Cf (JS, d ) x 
C^(M. dxd ). 

Proof. The proof follows exactly the steps of [3, Proposition 4.3], the only difference being 
that our Lemma |3~T1 has to be applied here in order to get the equivalent of inequality (28) 
in [3]. Details are left to the reader. 

□ 

With these preliminaries in hand, we can now turn to the finite dimensional distribution 
(f.d.d. in the sequel) convergence, which can be stated as: 

Proposition 2.5. Under the assumption \l.l[ let (X e ,X 2 ' e ) be the approximation process 
defined by (0|) and |7|). Then 

f.d.d. - lim(X £ , X 2 ' £ ) = (B, B 2 ), (9) 

where f.d.d. — lim stands for the convergence in law of the finite dimensional distributions. 
Otherwise stated, for any k > 1 and any family {sj, tf, % < k, < Sj < ti < T}, we have 

£ - hm( A7 x , X 2 ^ , . . . , Xt k , X 2 f tfc ) = (B tl , B 2 tl , . . . , B tk , B 2 ^. ) . (10) 

Proof. The structure of the proof follows again closely the steps of [31 Proposition 5.1], 
except that other kind of estimates will be needed in order to handle the Donsker case. 

To be more specific, it should be observed that the first series of simplifications in the 
proof of [3], Proposition 5.1] can be repeated here. They allow to pass from a convergence 
of double iterated integrals to the convergence of some Wiener type integrals with respect 
to X s . Namely, for % — 1, 2 and < u < t < 1, set 

Y\u,t)= ['(Bi-BDiv-uf-ldv, 

and for < u < t < 1 and (ui, . . . , n 6 ) in a neighborhood of in M 6 , set also 

Z u = Ul + u 2 B 2 u + u 3 Y 2 ( Uj t) + u 4 j (v - uf-^dW 2 v 

J u 



-U 5 



-u 6 



dw / (w-v) H ?({w-u) H 2 - (v - u)*) dW % 

J u 

dw / (w- v) H ~ 'Hw -u) H ^ dW 2 



Consider the analogous processes Y 1 ' 6 , Z e defined by the same formulae, except that they 
are based on the approximations 8 l,£ of white noise. We still need to recall a little more 
notation from j3j: for / G L 2 ([0, 1]) and t G [0, 1], we set 

d> £ (/) = E (e'JS/C^M**) , <jf f = ^ jf * f 2 (x)f 2 (y)I {lx _ yl<£ * } dxdy, (11) 



and 



$(/) = E (yfof^^i 



Then it is shown in [31 Proposition 5.1] that one is reduced to prove that \im e ^ Q v £ = 0, 
where v" is given by 



f £ = E ($> £ {Z £ )e 



£ \ P iw So ee,2 W rf « 



E (§(Z £ )e iw ti d£ ' 2< - u)du 



for an arbitrary real parameter w in a neighborhood of 0. Furthermore, bounding 
e iwf e e < (u)du trivially by 1 and conditioning, it is easily shown that v® is controlled by 
the difference E[\Q £ (Z £ )) — for which Lemma [3.31 provides the bound 

\§ £ {Z £ )) -HZ £ )\ 



< E 



./ |^ 3 (0^ £ )*||^ £ j| i2 A:3exp(4 W 2 ^||^||2 2 ) + w 2 ^ 2 - 1| ^ e j| a || ^ £ j| ^ 2 exp(w; 2 1| ^ e || | 2 : 



|2 > 
J e\\L 2 , 



+ - 7 = W 4 (0| e )^||Z £ ||i 2 ^ 2 exp(4 W 2 ^||Z e ||i 2 ) + - W 4 0| e exp( W 2 ||Z 

for any a G (0, 1). In order to reach our aim, it is thus sufficient to check the following 
inequalities: 

su P E(\\Z £ \\l) <M, lim£((0!J 2 )=O 
and for w < Wq, where Wq is a small enough constant, 

sup£(> 2|,Ze|1 ^ < M. 



However, these relations can be deduced, as (39), (40) and (41) in [3J, from Lemma [3.11 
(it should be noticed however that a one-parameter version of [21 Lemma 5.1] is needed 
for the adaptation of the latter result to our Donsker setting). The proof is thus finished 
once the lemmas below are proven. 

□ 



3. Moments estimates in the Donsker setting 

In order to deal with our technical estimates, let us first introduce a new notation: 
set pi = (1 — 5 x / 2 )/2 and p 2 = (1 + 5 1 / 2 )/2. Then the moments of any integral of a 
deterministic kernel / with respect to 6 t,£ can be bounded as follows: 

Lemma 3.1. Let m G N, / G L 2 ([0, 1]) 7 i G {1,2} and e > 0. Recall that the random 
variable rj is assumed to be almost surely bounded by a constant k v . Then we have 
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(12) 
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and 
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where M is the quantity defined at ( flip . 



Proof. We focus first on inequality ffT2l and divide this proof into several steps. 
Step 1: Identification of some key iterated integrals. Notice that 



I \ 2m 

/(r)^' £ (r)dr 



E 



< / |/(ri)| • • • \f(r 2m )\\E(e^( ri ) ■ ■■6^{r 2m ))\dr 1 ■ ■ ■ dr 2m . 
</[o,i] 2m 

Transforming the symmetric integral on [0, l] 2m into an integral on the simplex, and using 
expression (jl]) for 9 £ , we can write the latter expression as: 

(2m)! ^ ,w 2 



£ 2m 



E / ■•■/ -1/(^)1 



Kl , . . . , &2m — 1 
fcl>-"> fern 



X \ E (v l kl ■ ■■ik 2 J\hn>T2>->r2m} dT i ■ ■■dr 2m , (14) 
where n(e) = [^] + 1 and where we understand that f(x) = whenever x > 1. 

Let us study now the quantities E{rf ki - ■ -T] k ). If there exists / such that ^ ^ fcj 
for all j ^ / then E(rf ki - ■ -rjl) = 0. On the other hand, when k 2 i-± = k 2 i > k 2 i + \ 
for any /, we clearly have E[jf kx ■ ■ -rf k ) = 1. Finally, in the general case, for all / G N, 
l^ivli ' ' 'rfki)\ — K)- Separating the cases in this way for \E{rf ki ■ ■ ■^jt 2m )|, we end up with 
a decomposition of the form E[(J^ f{r)6 h£ {r)dr) 2m } = + T^, where 

(2m)! rkie 2 rkie 2 rk m e 2 rk m e 2 

— £ 2m E/ 

k\ , . . . , k m = 1 
k\ > ■ • • > fc m 



(fel-l)e 2 J(fei-1)£ 2 «/(fc m -l)£ 2 J(fe m -l)e- 



and where the term is defined by: 

(2m)!fc*» 



£ 2m 



• ■ ■ \f(r 2m )\ x /{ ri >r 2 >->r 2m }^i • ■■dr 2m , 



E ^i.-..,n., ( 15 ) 



rii, . . . ,n s > 2;s e {1, ...,m - 1} 
ni H + n s = 2m 



with 



= E / \f( r i)\---\f(r 2m )\I{r 1 >r 2 >->r 2m }dr 1 ---dr 2m , (16) 



ki , . . . , k s — 1 

&!>•••> fc s 



and where we have set D kl ... ka = n*=i[(^j — -O^ 2 ' ^i^ 2 ]™^ 

Let us observe at this point that we have split our sum into and because 
represents the dominant contribution to our moment estimate. This is simply due to the 
fact that is obtained by assuming some pairwise equalities among the random variables 



rf k , while is based on a higher number of constraints. In any case, both expressions 
will be analyzed through the introduction of some iterated integrals of the form 

If - 

K v (k;v,w) = — / \\f(r j )\I{ w > ri> ->r IJ >v}dr 1 - ■ -dr u , 

£ J l(k-l)e 2 M 2 }" j = i 

defined for u, k > 1 and < v < w < 1 . 

Step 2: Analysis of the integrals K v . Those iterated integrals are treated in a slightly 
different way according to the parity of v. Indeed, for v = 2n, thanks to the elementary 
inequality 2ab < a 2 + b 2 , we obtain a bound of the form: 

n(s) 

v,w) (17) 



k=l 

sr f ^ ( f 2 (^i) + f 2 (x 2i+1 ) \ T 

h ^ W,^- 1\ \ 2 ) w>->«->«}^i • • • ^ 

yj-JT / / 2 («l) • • •/ 2 (^n)^{ lU > :r i>-> a; „>j;}^l • • -^2n 

fc=l 5 -/[(fc-l)£ 2 ,fce 2 ] 2 ' 1 

"•(e) 

= f( x i) ■ ■ ■ f{x n )I {w > Xl >...> Xn > v} dxi ■ ■ ■ dx n 

J[(k-l)e 2 ,ke 2 ] n 

< / / 2 (£i) ■ • ■ f 2 (x n )I^ xl _ Xn<£ 2 } I {w > Xl ^...> Xn > v }dx 1 ■ ■ ■ dx n . 
J\o,i\ n 



< 

n(e) 



The case v = 2n + 1 can be treated along the same lines, except for the fact that one has 
to cope with some expressions of the form 

n(e) n(e) 



y~]K 3 (k;v,w) < y~]- I \f(xi)\f 2 (x 2 )I{ w > Xl > X2 > v }dx 1 dx 2 

k =i k =i e J[{k-iy 2 M 2 ? 

< - / \f(x 1 )\f 2 (x 2 )I{ Xl - X2< e2}I{w>x 1 >X2>v}dx 1 dx 2 . (18) 

£ ^[o,i] 2 

Combining ( {TBI and ( fTTj) we can state the following general formula: let v > 1, and 
define a couple (z/*, P) as: (i) z/* = u/2, v = if v is even, (ii) z/* = (r/+ l)/2, z> = 1 if v 
is odd. With this notation in hand, we have: 

"(e) p 

J2K„(k-v,w) < - / \f(x 1 )\ 2 -»f(x 2 )---f 2 (x„*)I {xi ^ <£ > } 

x /{«,>xi>->x„.>:«}rfa;i • • -c/av*. (19) 

S^ep 5 1 .- Bound on T^. It is readily checked that can be decomposed into blocks of the 
form K 2 (k;w,v), for which one can apply ( TL9l ). This yields 

(2m)' n(£) /" fcl£2 f km£2 



Or 

< ( 2m ) ! nn|2™ 



Step 4 : Bound on U ni na . Recall that U ni na is defined by (1161) . We introduce now a 
recursion procedure in order to control this term. Namely, integrating with respect to the 
last n s variables, one obtains that 

^ j 2m-n s 

-^U ni ,...,n a = £ 2m-n s E /„ II l/( r U{ri>r 2 >--->r 2m _„ s } 



fci, . . . , fc s _i = 1 
fci > • • • > fc s _i 



^fcl-fcs-l z=i 



x K ns {k s ;0; r 2m. ) rfri • • • dr 2m 

— n s 

Plugging our bound (JTHJ) on i^„ s into this expression, we get 

1 -, "(e) 2m— n a 

^ U n^n s <-^~ £ / / . II l/Cn)| l/Cl/!)! 2 -*- 

c , JD k ,,.. k _ , j 0,1 n « i_i 



fcii ... , fc s _i = 1 
ki> ■■■> fc s -i 



x n/ 2 (%)^i-y n| <e 2 } J {ri>r 2 >--->r 2m _ ns > yi > - >^ ■ • • dr 2m - na dy 1 ■ ■ ■ dy n * . 
i=2 

We can now proceed, and integrate with respect to the variables r\ for < I < k s . In 
the end, since n s = 2m, the remaining singularity in e is of the form Y[ £~ ns ■ However, 
each of the singularity e~ na comes with an integral <f> e with q = If] 1 ^ 2 . The latter integral 
is easily seen to be of order e, which compensates the singularity e~ ns (recall that h s < 1). 
Hence, iterating the integrations with respect to the variables r, we end up with a bound 
of the form 



1 



^u ni ,..., ns < ^^y\\f\\T- 2) ff < (^rll/ll^ra 1 - ( 20 ) 

Step 5: Bound on T^. Owing to inequality f[2"Ul) . our bound on can be reduced now 
to an estimate of the number of terms in the sum over ni, . . . , n s in formula (|15p . This 
boils down to the following question: given a natural number n, how can we write it as a 
sum of natural numbers (larger than one)? 

This is arguably a classical problem, and in order to recall its answer, let us take a 
simple example: for n = 6, the possible decompositions can be written as {6; 2 + 2 + 2; 2 + 
4; 4+2; 3+3}. Furthermore, notice that the decompositions of 6 can be obtained by adding 
+2 to the decompositions of 4 or adding 1 to the last number of the decompositions of 5. 
Extrapolating to a general integer n, it is easily seen that the number of decompositions 
can be expressed where (u n ) n >\ stands for the Fibonacci sequence. We have thus 

found a number of decompositions of the form 

N n = s- 1 / 2 1 - P r n , 

where the quantities pi,p2 appear in formula (TT2T) . Moreover, the number of terms in T 2 
is given by A^ m — 1> the —1 part corresponding to the term 7\. 

Putting together this expression with ( [20]) and the result of Step 3, our claim ( fT2l) is 
now easily obtained. 

Step 6: Proof of ( T7^) . The proof of ( [TBI follows the same arguments as for ( fT2i) . We 
briefly sketch the main difference between these two proofs, lying in the analysis of the 



term U nit4 __ >ns . Indeed, since we are now dealing with an odd power 2m + 1, the equivalent 
of ( l20|) is an upper bound of the form 



{\yi-y2\<e 2 }dyidy2 \f { x j)I{xi>—>x m -i}d x i """ dx m —i. 

£ Jo Jo Jfo.n™- 1 TJT 



[o,i] r - - j=1 

Furthermore, applying Holder's inequality twice, we obtain 



(21) 



1 r i 



JO 



(y2)I{\ yi -y 2 \<e 2 }dyidy 2 = - f (y 



!A(y 2 +£ 2 ) 



0V(j/ 2 -£ 2 ) 



\f(y 1 )\dy 1 dy 2 



< 



1 r l 



JO 



/ (yi)f (y 2 )I{\ yi -y 2 \<e*}dyidy< 



l 2 = WfrMUv 



and thus we can bound f[2"Tj) by ( - m i 1 - )! (4>f) 2 11/11 which ends the proof 



( m _l)i iv/ 



□ 



Our next technical lemma compares the moments of a Wiener type integral with respect 
to 9 £ and with respect to the white noise. 

Lemma 3.2. Let m G N, / G C a ([0, 1]), i G {1, 2}, e > and for m > 1, sei 

2m" 



(2m) 



/(r)0 iie (r)dr 



2 m m\ 



f(si) ■ ■ -/ 2 (s m )rfsi •• - ds r 



[0,1]' 



T/ien 

We W Ji<£ 2a ||/|U 2 ||/|| Q . 

^ -For any m > 1, £/ie following inequality holds true, where we recall that pi,p 2 have 
been defined just before Lemma \3.1[ 



1 



_2a 



fm - 1)! 



\2m-l 



a\\J II L 2 



1,2m 



2m-l J2m-1\ ( ie\~ II f\\2m-2 



1 



v/5(m-2)! vri r2 " JIIL/ (m-2)\ 

Proof. We divide again this proof into several steps. 

Step 1: Variance estimates. We prove here the first of our assertions: Notice that 

n ( £ ) r ke 2 rke 2 



2(m-2) , £ 

L 2 



? / / 2 (^i = ^E / / / 2 (*i)^ 1 . 

2 JO A£ k=l J(k-l)e 2 J(k-l)e 2 



On the other hand 



2 



1 \ 2' 

f(r)6^(r)dr 



^ /-fce 2 /.fee 2 

I I f(n)f(r 2 )dr 2 dn. 



k=l 



(fc-l)e 2 J(A;-l)e 



We thus get 



Ji 



-E/ / f(n)(f(r 2 )-f( ri ))dr 2 dr, 



2e 2 



k=l 



(fc-l)e 2 J(k-l)e 



2i2 / / o /(ri)(/(r 2 )-/(rx)) EwgD ^ 



and hence this quantity can be bounded as follows: 



Ji < 



< 



2e 2 
1 



1 rl 



JO 
1 



|/(ri)||(/(r 2 ) - /(^l))|/{|r- 1 -r 2 |< £ 2}^2<iri 
a I ^1 1 I{\ri— r2\<e 2 }d'T2dTi 



< e 



2a 



L' 2 \\J a 



which is the first claim of our lemma. 

Step 2: decomposition for higher moments: We can follow exactly the computations of 
Lemma l3~T| Step 1, in order to get 



1 



(2m)! 



E 



l \ 2m 

/(r)^' £ (r)dr 



o 



m ' mi 



w * h f L = (2% 
and we obtain 



for j = 1, 2. Furthermore, the term can be bounded as in Lemma 1371 



irii < 



1,2m 



V5(m - 2) 



2m-2 
L 2 ■ 



(22) 



Step 3: Study ofT^: We analyze in a slightly different way as in Lemma I3TT1 Namely, 
we first write 



7; 1 



"i 2 m e 2m 



n(e) 

£ 

k\ , . . . , k rn — 1 

k\ !> • • • !> fc m 



(fci-l)e 2 J(fci-l)e 2 ^(fc m -l)e 2 ^(fe m -l) £ 



/(n) ■■■f(r 2m ) 



^{{n,r 2 }>->{r2m-i,r 2m }}dri • ■ ■ dr 2m , (23) 



where we have written {a, 6} > {c, d} for a A 6 > c\/d. We will now compare this quantity 
with another expression of the same type, called and defined by 



fi 



2 m m! J m 



f{si) - ■ ■ f 2 {s m )dsi - ■ - ds f 



Let us thus write Ti, as 



If m 

Z J[0,l] m j=1 



J(fci-l) 6 a </(fc m -l) £ 2 

M i ■ ■ ■ i — 1 

^1 ^ * * ' ^ km 

n ( £ ) /-fcie 2 /-fc m£ 2 m 



^ / "/ T/ 2 ( S j) I {si>->Sm} ds l--- ds m + Tl l , 



2 

fel , . . . , fc m = 1 

fci > • • • > k 



where represents the part of the sum taken over the indices k\, . . . , k m such that there 
exist I satisfying k\ = ki + \. However, this latter term can be bounded as in ( 1201) . yielding 

Step 4 : Conclusion. Putting together the decompositions we have obtained so far, we end 
up with 



Jm < \fl\ + |T£ | + 



n ( £ ) r k l£ 2 r k ie 2 r k m e 2 r k m e 2 
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Invoking our estimates ( l22l) and ( l24l) on and T^, our bound on J m easily reduced to 
check that 



1 



2 m m\e 2m 



|/(ri)---/(r ara )-/ 2 (ri)/ a (rs)---/ a (ra m -i)| 



[0,1] 



2q|| rll II r\\2m-l 
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The latter inequality can now be obtained from the decomposition 

\f(r 1 )---f(r 2m )-f(r 1 )f(r 3 )...f(r 2m . 1 )\ 
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the inequalities 
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L 2 ' 
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and from the estimate we have already obtained for J\. This finishes the proof. 



□ 



Finally, the characteristic function of a Wiener type integral of the form f(r)9 k,£ (r)dr 
can be compared to its expected limit f f(r)dW k in the following way: 

Lemma 3.3. Let f e C a ([0, 1]) for a certain a 6 (0, 1), k 6 {1, . . . , d} and e > 0. For 

any u G R, we have: 



E [ e iu fo /(r)* fc ' e (r)drj _ £ fgiu J,, 1 /(r)d^ j 



< 4(l/5) 1 /V(0 e f )^||/|| L2 ^exp(4 M 2 ^||/||i 2 ) +u 2 e 



L 2 



exp(w 



2|| ^||2 > 

L 2 1 



Proof. Let us control first the imaginary part of the difference. Using lemma 13.11 and 
invoking the fact that the odd moments of a Gaussian random variable are null, we get 



Im ( E [e iu f(r)o h ' e (r)dr'j _ E ft f(r)dw r ^ 
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m=l ^ 

<4(l/5) 1 /V(0})^||/|U 2 ^exp(4 M 2 ^||/||i 2 ). 
In order to control the real part of the difference, we will use Lemma 13. 21 This yields: 
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The latter quantity can be bounded by 



u 




which ends the proof. 

□ 
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